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We consider magnetic oscillations of resistivity of a clean (mean free path l R) hollow cylinder 
with fluctuating (with an amplitude of fluctuations A R <C R) radius R, threaded by magnetic flux 
d>. We demonstrate, that for weak fluctuations (A R <C pf 1 ) the oscillations have a standard period 
2$o, characteristic for oscillations in a clean system, while for A R 3> pf 1 they become 'f’o-periodic, 
which was expected only for dirty systems with l -C R. The work is motivated by observation of 
predominantly ^o-periodic magnetic oscillations in very clean Bismuth wires. 


Introduction. Oscillations of resistivity p of a hol¬ 
low conducting cylinder with change of magnetic field 
H , parallel to the axis of the cylinder, is a well known 
phenomenon, manifesting the quantum nature of con¬ 
ducting electrons in metals. The basic reason for these 
oscillations is the interference of electrons traversing the 
cylinder along paths with different topological winding 
numbers. In the simplest case of clean cylinder (with the 
circumference 2nR <C l) it is actually a direct realization 
of the Aharonov-Bohm effect |lj ; the resistance is a peri¬ 
odic function of the flux $ = t:R 2 H with the period 2$o, 
where $o = 7rc/i/2e is the flux quantum (see, e.g., 3. 4]). 
The general arguments (c.f. i) show that the same pe¬ 
riodicity should be observed for any thermodynamic or 
kinetic property of a cylinder - for arbitrary relation be¬ 
tween 2nR and l. In particular, for the resistivity 
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In dirty cylinders, where 2nR l the odd-n harmon¬ 
ics of the AB-oscillations are washed out due to strong 
variations in the length of different diffusive trajectories 
that lead to randomisation of the non-magnetic part of 
the phases of electronic wave-functions. In a seminal pa¬ 
per by Altshuler, Aronov and Spivak [5] it was shown, 
however, that the even-n harmonics - the oscillations, 
associated with a special sort of trajectories (the ones, 
containing closed topologically non-trivial loops on the 
cylinder) survive the randomisation. These even-n oscil¬ 
lations arise due to the weak localisation corrections and 
manifest the interference between the contributions of 
paths, corresponding to the same trajectories, but with 
opposite directions in which the loop is traversed. In 
such an interference the randomised nonmagnetic part 
of the phase does not enter, while the magnetic part of 
the phase is doubled and, therefore, the AAS-oscillations 
occur with the period $ 0 . The AAS-oscillations were 
observed experimentally (see [§, ]}), although their am¬ 
plitude is small: ~ (ppl) -1 - 

Thus, a simple correspondence seemed to be estab¬ 
lished: In clean cylinders with 2ttR <C l the AB- 
oscillations with the period 2<I>o should dominate, while 


dirty cylinders (with 2ttR l) is the domain of the AAS- 
oscillations with the period $ 0 - 

Recently, however, the dominance of the <f>o-periodic 
(AAS-like) oscillations was experimentally observed Q 
in a manifestly clean Bismuth wires, where the AB- 
oscillations would be natural to expect. Note, that in 
this system the surface states are believed to be mainly 
responsible for the conductivity and the spin-orbit inter¬ 
action is important for these states. In this situation the 
oscillations of resistivity were discussed in [§j, where it 
was shown that in the case of strong surface disorder the 
odd harmonics in m should be suppressed, and the os¬ 
cillations should be indeed $ 0 -periodic, like in the AAS- 
case. However, the model, used in @ seems to be not 
applicable to the experiments Q, since it implies strong 
disorder, which is apparently absent there. 

Trying to explain thie unexpected experimental obser¬ 
vation of [§j, one should address the following question: 
Suppose that there is some mechanism of suppression of 
odd-?i harmonics, then how is it possible, that the same 
mechanism does not suppress also the mean free path l 
and does not show up in the monotonic part of the resis¬ 
tivity? 

On the other hand, under the condition 2nR <C l 
the standard method of finding even-n AAS-oscillations, 
based on the diffusion approximation, does not apply, 
and even the qualitative physics underlying the 'bo- 
periodic oscillations may be different from the weak lo¬ 
calisation effect, considered in [jjj]. 

In this Letter we propose a simple mechanism for this 
effect, related to weak and smooth fluctuations of the 
radius of the wire. Here we will discuss it only for the 
simplest case of weak short range impurities, that can be 
treated within the Born approximation; we also do not 
take into account the spin-orbital interaction. The latter 
effects will be considered in a separate publication, here 
we only mention, that they may, in some cases, consid¬ 
erably alter the amplitude of the even-n harmonics, but 
they do not change the principal result: odd harmonics 
are suppressed much stronger than even ones. 

We argue, that even relatively small fluctuations of 
the cylinders radius R{z), measured at different cross- 
sections z, are able to exponentially suppress the oscilla- 


2 


tions. Namely 
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A n (Q, q) = a n -\ - A_ = eT 2n ' 1 ^ cos(2 nQ + 7r/4) 

V7rnQ 

_ 6 f 0, (ra - odd), 

" tt (nQ) 3 ^ 1, (n — even), 

where we have introduced the parameters 

Q = nR pf, q = nARpp <C Q, (3) 

R being the cylinders radius, averaged over the positions 

/-=—\ 1/2 _ 

of the cross section, and A R = I (R— R) 2 J -C R, be¬ 
ing the relatively small variance of R. The quasiclassical 
parameter Q we will assume to be large: Q 1, while q 
may be either small, or large. 

We demonstrate that at q 3> 1 the fluctuations expo¬ 
nentially suppress all the amplitudes A n ; however, while 
the odd harmonics are suppressed to zero, the even ones 
saturate at some small (in parameter Q ), but finite 
values a n . Thus, at q —> oo the oscillations become 
periodic, as in the dirty case. 

On the other hand, if the fluctuations of R(z) are adi- 
abatically smooth (their correlation radius £ 3> R) then 
this kind of disorder leads to only exponentially small 
contribution to the scattering of electrons and, therefore, 
practically does not affect the mean free path l. 

It should also be noted that our result (0, in contrast 
with the result of Q , is not proportional to small param¬ 
eter (ppl)~ 1 '- indeed it does not involve any weak locali¬ 
sation effects and does not depend on l at all. Although 
the derivation of expression for the residual amplitude a n 
requires taking into account quantum corrections to the 
standard Drude formula, these corrections are not asso¬ 
ciated with the interference of different scattering events. 
It does not mean that the localisation corrections do not 
exist in our case - they do, but for low concentration of 
impurities the corresponding corrections to A n are small. 

A homogeneous tube. We start with considering a 
hollow metallic cylinder with a fixed radius R and length 
L, threaded by a magnetic flux $ (the magnetic field is 
oriented along the axis of a cylinder z ). The effective 
two-dimensional metal, constituting the tube, is charac¬ 
terized by an isotropic spectrum of electrons, pf and vf 
being, respectively, the momentum and the velocity at 
the Fermi level Ep- There are rare impurities in the 
system - with the two-dimensional concentration nj nip , 
the corresponding mean free path being l. In this paper 
we restrict our consideration to the case of low concen¬ 
tration n^p < R~ 2 and consider the semiclassical and 
clean case, when 

Pp 1 < 2? tR < Azi rnp , l <C L (4) 

where Azi mp ~ {2irRn^ p )~ 1 R is a characteristic 
separation between impurities. 


In an absolutely clean system the eigenfunctions of the 
hamiltonian at the Fermi surface are 

z) = exp {ik m z + imcj)}, (5) 

k m = ± \JPp — (m + $/2$o ) 2 /R 2 , (6) 

where m € Z is the azymuthal quantum number. 

Under the condition 0 the system can be treated as a 
multichannel quasi-one-dimensional one (the number of 
channels N c h = (2ttR)(pf/ 2n) = pfR 1), so that the 
localization length £i oc = N^l/2 3> l is large. We will 
assume that it is larger than the length of the wire: 

Aoc = HjpfR)/ 2 (7) 

and the specific one-dimensional localization effects are 
relatively small, and can be treated perturbatively (weak 
localisation). From this assumption it immediately fol¬ 
lows, that, in the leading semiclassical approximation, 
the conductivity of the wire a^ ld ^ can be described by 
the standard Drude formula 

= {2TrR^° ) )e 2 D 0 = e 2 C loc , D 0 = v F l/ 2, (8) 

where Dq is the diffusion coefficient, and = pf/2ttvf 
is the density of states of an infinite two-dimensional 
metal at the Fermi surface (without “spin two”). The 
factor 2pRvp^ = ( pfR)/vf = N^/vf is nothing else 
but the effective one-dinrensional density of states in the 
wire. 

Certainly, there is no place for the dependence on $ in 
the result ©, such a dependence can only show up in the 
corrections to the semiclassical approximation. There are 
three types of these corrections: 

1. “Thermodynamic” corrections. Here we gather all 
the oscillatory effects, that arise already in the ideal tube 
without any impurities at all - corrections to the density 
of states and similar. In the leading order in the semi¬ 
classical parameter Q > 1 the largest of these corrections 
are ~ Q~ 1 7 2 <C 1, but we will need also the smaller ones. 

2. “Kinetic” corrections - they arise due to oscilla¬ 
tory corrections to the single-impurity scattering ampli¬ 
tude. Usually the single-impurity effects only lead to 
an irrelevant renormalization of parameters of an impu¬ 
rity. In our case, however, due to nontrivial topology 
of the system, these corrections are highly sensitive to 
the magnetic flux 4> and exhibit oscillatory behaviour. 
The relative amplitude of kinetic corrections (as well as 
the thermodynamic ones) does not depend on the con¬ 
centration of impurities, its smallness comes mainly from 
the semiclassical parameter. The kinetic corrections are 
proportional to the scattering amplitude, they arise only 
beyond the Born approximation. In this Letter we will 
restrict our consideration to the first Born approxima¬ 
tion, so that the kinetic corrections will not be taken 
into account in what follows. 
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3. Localization corrections. They stem from the effect 
of interference of scattering events at different impuri- 
ties, their relative amplitude is proportional to n im and 
is small in parameter l/(ppl) -C 1. The localisation cor¬ 
rections will also be neglected in this Letter. 

Conductivity in the Born approximation. Con¬ 
sider randomly placed short range impurities with con- 

/Q\ 

centration (that is, the average number of impurities 
per unit area of the cylinders surface). Their potential is 

tfimp = V6(r - r i) (9) 

i 

where ri (z = z*, (f> = <pi) is a position of i-th impurity 
on the cylinder. In this paper we consider only the case 
of weak impurities, with small dimensionless strength 


where v m k = v F k/pF is the velocity in the z direction 
and A is, again, the correction, arising due to the finite 
radius of the cylinder. 

Thermodynamic corrections. Corrections A and 
A are the only sources for the oscillations of the resis¬ 
tivity within the Born approximation. Here we calculate 
them with the help of the Poisson summation formula. 

An exact expression for vp reads 

v F = 5 ( e f - [k 2 + (to + <f>/2<f> 0 ) 2 /R 2 ]/2M) = 

mk 

oo C rft 1 

= E 

n=— oo J ' 

x exp {2nin(pR sin 9 + 4>/2$o)} , 

(16) 


A = 7Ti4 0) V<l (10) 

so that the Born approximation is valid. 

For point impurities the scattering is isotropic and 
therefore the transport time coincides with the simple 
decay time. The latter for a state (mk) can be written 
with the help of the Fermi golden rule 

r Km = 27m imp E \ V kk'mm'\ 2 S[E(kni) - E(k'm')\ (11) 

k'm' 

Since for the potential © matrix elements \Vkk'mm'\ 2 
depends neither on km, nor on k'm ', we can write 

T k,m = 2 ™imp (12) 


where 


f f = Y, S ( Emk - = v f\ 1 + A )> ( 13 ) 

mk 

is the density of states at the Fermi level for an ideal 
tube. Here A is the correction, arising due to the finite 
radius of the cylinder. 

It is important to note, that the decay rate in the lead¬ 
ing approximation (i.e., without taking into account the 
correction A) is isotropic: it does not depend on the an¬ 
gle 9 between the momentum and the cylinder axis. 

Within the Born approximation one can write for the 
conductivity 

<7 = e 2 (2pR) v 2 (mk)Tt r 5[E(km) — Ef] = 

mk 

= e 2 ^_ F( 27r r ~)~ f ( 14 ) 

where the “transport density of states” is defined as 
ff = 2j2(k/p F ) 2 S(E mk - E f ) = 4° } (1 + A) (15) 

mk 


where we have introduced effective mass M = pf/vf 
and two-dimensional momentum: p z = k = pcosO and 
Px = (m+<f>/2<l> 0 )/R = psind, 9 being the angle between 
the momentum and the cylinders’ axis. Performing a 
trivial integration over p and introducing parameter Q = 
ttpfR 1, we get 


{ °°^ r'n 1 

1 + 2 V / -—cos(2?r(5sind) cos(7m$/$o) ) 

7^i 2?r J 


Thus 


OO 

A = 2 ^ J 0 (2nQ) cos(7m$/$ 0 ) (17) 

n =1 


where Jo is the Bessel function. Using the two leading 
terms of its asymptotics, we get 


OO 


A = ^ {A(2nQ) cos(2nQ — 7r/4)-f 

n=1 

+ B(2nQ) sin(2n<5 — 7r/4)} cos(7m<I>/<f>o) 


(18) 


where 

^> = 2 \/5- SW = 2 \/E (19) 

In the leading approximation in Q^ 1 ! 2 


A = E 


2 

(n Qn) 1 / 2 


cos(2 Qn — 7t/ 4) cos(7m$/$ 0 ) 


( 20 ) 


but we will also need the subleading corrections in what 
follows. 

For the transport density of states we have 


v F = 2 ^2(k/p F ) 2 

mk 

xd (E f - [k 2 + (to + $/2$o) 2 /i? 2 ]/2M) = 

= 4 0) (1 + A), (21) 
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where 


OO 


A = 2 YJ, {</o(2nQ) + J2(2,nQ)} cos(7m$/$ 0 ) 

71=1 

• (22) 

Using the asymptotics of the Bessel functions, 
finally write 

we can 

OO 

A = 'Y / jA(2nQ) cos(2nQ — 7r/4)+ 

71= 1 


+B(2nQ) sin(2n<5 — 7r/4)| cos(7m4>/<h 0 ); 

(23) 

nr 3 ~ nr 2 

A(x) = 2J — —, B(x) = 2J - 

V 7ra; 4ar \ nx x 

(24) 


Thus, we conclude that the oscillation correction to vf is 
much smaller, than that for vp\ A<C A. The reason is 
that the oscillations are provided by the states with very 
small jfe’s, whose velocity along z is very small. 

Now we are prepared to write a expression for the re¬ 
sistivity 


1 

P 


l + A 

1 +A 


e 2 D 0 (27rR)F { °\ 


(25) 


2A 2 n 


( 2 ) 


,(°) 


D 0 = 


TpVp 
2 ’ 


(26) 


are the perturbative decay rate and the diffusion coeffi¬ 
cient of an infinite two-dimensional metal. 

Corrections to the resistivity. The relative oscil¬ 
lations of the resistivity in the second order in A are 


A 


(i) 

tot 


6p A-A (i) 

P 1 + A tot 

_ a — A A*- 2 ) - 


+ A (2 ^ -)-••• 

I ^tot V , 

= -(A -A)A 


(27) 

(28) 


In the linear term the contribution A always dominates, 
so that A« « A is given by formula ©. 

/o') 

In the second order correction A( 0 ( we should 
keep the “diagonal” terms of the form cos 2 (2?rQ — 
7t/4) cos 2 n(7T?r$/$ 0 ) or sm 2 (2nQ — 7r/4) cos 2 (7m$/<l>o)- 
Only these terms contain contributions (equal to 
(1/4) cos(27rn4j/<ho)), that are insensitive to the fluctua¬ 
tions of the cylinder radius R , since they do not contain 
the factors, oscillating in Q: 


A^ = - (A - A)A « F(2nQ) cos(2t rn$/$ 0 ) (29) 

n= 1 


F(x) = ~{{A-A)A + (B-B)b} = -^ (30) 

In a homogeneous tube the correction A^ should be 
neglected, compared to A^, but, as we will see, it can 


become dominant, if fluctuations of radius are taken into 
account. 

The series (l29l) converges for all $, therefore the cor¬ 
rection A£ 2 j is always relatively small due to the factor 
Q~ 3 . It is not the case of the correction A^, since the se¬ 
ries m demonstrates a square-root divergence at certain 
values of <f>. Below we discuss the behaviour of the resis¬ 
tivity close to the singularities and explain their physical 
origin. 

Singularities in p(<E>) . It is convenient to rewrite the 
result CU1) in the form 

A = (ttQ)- 1 / 2 

xRe {S(2nQ + 7t4>/$ 0 ) + S(2nQ - tt$/$ 0 )} (31) 

where 


*-50 A Ti 7 

(32) 

i V n 

n= 1 

This function has a singularity at z —> 2nN for any 
integer N. In the vicinity of the singularity, where 
8z = z — 2 ttN < 1, one can write 


S{z) 


= e“ i7r/4 Y 

n =1 


ginSz 




dn AuSz _ 


g-m/4 r oo e -t dt _ r^ 

\J—iSz Jo y/t \ Sz 


(33) 


Thus, S(z) is real and positive for <5^: > 0 and it is purely 
imaginary for Sz < 0. Then 


A » [27rQ(Q/7r ± $/2$ 0 - A)]” 1/2 

x0(Q/7t±$/2$o- N) (34) 


The physical reason for this square-root divergence is as 
follows. In our ideal quasi-one dimensional system we 
have a series of one-dimensional bands and, when the 
Fermi energy Ep crosses the edge Epf of one of these 
bands, the density of states diverges (as it should be in 
one dimension!): 

v F = Y V N d \ E F )> 

N 

( (0) \ 1 / 2 

y ( N d) (Ep ) ~ f e ” f _ Ef j 0{E N - Ep) (35) 


A »- =d(E N - E f ) (36) 

2fR\J 2iri'p' > (En — Ep) 

A natural question arises about a possible mechanism for 
smoothing the singularities in the p( ( I , )-dependence. The 
most obvious one is smearing of the singularities already 
in the density of states vf. Such a smearing should occur, 
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when the wave-length of an electron in a corresponding 
one-dimensional band becomes of order of the mean free 
path l. At low concentration of impurities, however, l is 
very large and another mechanism of smoothing seems to 
be more effective - the non-Born corrections to the scat¬ 
tering amplitude, which dramatically increase close to the 
singularities and lead to effective suppression of scatter¬ 
ing probability, inspite of the divergence of the density 
of final states. This is a completely single-impurity ef¬ 
fect, so it does not depend on the impurity concentration 
(and, therefore, on l ) at all. 

All this singularity story is only actual for a homoge¬ 
neous tube, however. As we will soon see, the fluctua¬ 
tions of the radius destroy the singularities anyway, and, 
therefore, in this Letter we will not discuss the singularity 
physics in detail and postpone it for a separate paper. 

Fluctuations of the radius. Now we consider a hol¬ 
low metallic cylinder with fluctuating radius R(z). For 
simplicity we will assume that the cylindrical symmetry 
is conserved, though this assumption is not essential: it 
is only the local circumference c(z), that matters, so that 
one can define effective local radius, as R(z ) = c(z)/ 2n. 
The fluctuations of R{z) are taken in a gaussian form: 

V{R(z)} = exp | — - J dzdz'Q(z — z')8R(z)5R(z')^ , 

8R(z) = R{z) - R, (SR{z)SR(z')) = K(z - z'), 

I< = k- 1 = (A Rff(\z - *'|/0, 

f(x) being some function with /(0) = 1 and characteris¬ 
tic scale ~ 1. The fluctuations are supposed to be rela¬ 
tively small and smooth: 

AR < R, £ > R, (37) 

So smooth fluctuations only lead to an exponentially 
(in parameter £kp) small scattering of electrons and thus, 
do not give any considerable contribution to the decay 
rate. On the other hand, as we will see, they strongly 
suppress the Aharonov-Bohm oscillations. 

For simplicity we discuss only the case of very long 
correlation length: l -C £ <C L. In this case one can 
think about our system as a sequence of short subwires 
with lengths Li < £, each of them having its own radius 
Ri, randomly distributed according to the Gaussian law. 
Each subwire has its own resistance IZj = p(Rt)Li and 
the total resistance 1Z of the system of subwires, con¬ 
nected in series, is then 

n = 53 = E P( R i)k = L k (38) 

~ p = J 7^ ap { J ~^§~} p{R>dR (39) 

Obviously, the averaging will dramatically suppress the 
Aharonov-Bohm oscillations, if the amplitude of fluctu¬ 
ations A R Pp 1 . Substituting (1271) (1751) and (1701) to 


(1391) . and performing trivial integration, we arrive at the 
final result ©■ 

Conclusion. We have shown that the character of 
oscillations in the p(<l>) dependence for a cylinder with 
the fluctuating radius is controlled by the parameter q = 
7tA Rkp, where AR is the amplitude of fluctuations. 

1. For q < Cl the oscillations are 2$ 0 -periodic, highly 
unharmonic, with sharp asymmetric maxima: 

/ 00 2 — 

P = Po 1 + E _ cos(2nQ + 7 t/4 ) cos(7r?r$/$ 0 ) 

\ n= 1 V' Kn Q 


The relative amplitude of oscillations at maxima depends 
on details of the impurity potential and may be large. 

2. For q ~ 1 the oscillations remain 2<f>o-periodic, but 
become almost harmonic and have small amplitude A\ ~ 


2 _ 9 2 

1 -I-^=e ~ q cos 

V7 tQ 

3. For g > 1 the oscillations become c E > o-periodic, al¬ 
most harmonic, and have still smaller amplitude Ai ~ 

Q~ 3 < A-r- 


(2 Q + 7t/4) cos(7t4>/$o) 



p ~ Po ( H- E cos(2tt$/4> 0 ) 

V 2t tQ 

The $o-periodic oscillations are due to interference of 
oscillating corrections to two different thermodynamic 
quantities: to the standard density of states and to the 
“transport density of states”. Therefore the effect has 
additional semiclassical smallness (higher power of Q ). 
It is important to stress, however, that, in contrast with 
the Altshuler, Aronov and Spivak effect, our one does not 
involve the interference between the scattering events at 
different impurities and, therefore, its relative amplitude 
does not contain small parameter (pfI ) 1 and does not 
depend on the concentration of impurities. 

The author is indebted to M. V. Feigel’man, P. A. Iose- 
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